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Abstract 

The middle levels conjecture asserts that there is a Hamiltonian cycle in the middle two 
levels of 2k + 1-dimensional hypercube. The conjecture is known to be true for k < 17 
[I. Shields, B.J. Shields and CD. Savage, Disc. Math., 309, 5271-5277 (2009)]. In this 
note, we verify that the conjecture is also true for k = 18 by constructing a Hamiltonian 
cycle in the middle two levels of 37-dimensional hypercube with the aid of the computer. 
We achieve this by introducing a new decomposition technique and an efficient algorithm 
for ordering the Narayana objects. 

1 Introduction 

Let Q n denote the n-dimensional hypercube, i.e., Q n is a graph with 2™ vertices, 
each vertex is labeled by an n-bit binary string and two vertices are adjacent iff 
their strings differ exactly in one bit. The i-th level of Q n is the set of vertices 
labeled by strings with exactly i ones. 

The middle levels graph is a subgraph of Q2fc+i induced by the middle two levels 
k and k + 1, and is denoted by M 2 k+i (see Fig. [[]). The middle levels conjecture 
asserts that the graph M2k+i has a Hamiltonian cycle for every k. It appears as an 
"exercise" in Knuth's book [3l Exercise 56, Sect. 7.2.1.3], in which the conjecture 
is credited to Buck and Wiedermann pQ. 

In spite of considerable efforts, the conjecture remains open (see e.g., [2j d] and 
the references therein). It was shown to be true for k < 11 by Moews and Reid, 
and for 12 < k < 15 by Shields and Savage [5] and 16 < k < 17 by Shields et al. 

In this note, we verify that the conjecture is also true for k = 18 by constructing 
a Hamiltonian cycle in the middle two levels of 37-dimensional hypercube with the 
aid of the computer. We achieve this by plugging a new decomposition technique 
and an efficient algorithm for ordering the Narayana objects into a Hamiltonian 
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path heuristic developed by Shields et al. [UH]. In the largest case, our program 
could find a Hamiltonian path in a graph with ~ 2.36 • 10 8 vertices in about a week 
on a standard PC. 

The organization of this note is as follows. In Section [21 we briefly review the 
approach taken by Shields et al. [U [5] for reducing the size of the problem. In 
Section El we describe an additional reduction that decomposes the problem into a 
number of smaller subproblems. In Section^ we introduce an efficient algorithm for 
ordering the Narayana objects which was helpful for reducing the resource needed 
in the computation. Finally in Section 5, we summarize our computational results. 
Throughout the paper, n = 2k + 1 denotes the dimension of a hypercube. 

2 Reducing the problem 

The Hamiltonicity of the middle levels graph, which has 2Q) vertices, can be 
reduced to the problem for finding a suitable Hamiltonian path in a smaller graph 
with (2) I n vertices [5] . 

For an ra-bit binary sequence x = x\X2 ■ ■ ■ x n , define the cyclic shift a by a(x) = 
X2X3 ■ ■ ■ x n x\. For every two vertices x and y in M n , x and y are adjacent iff a(x) 
and a{y) are adjacent. This naturally introduces an equivalence relation ~ on the 
set of vertices of M n such that x ~ y iff x = a l {x) for some integer i. By noticing 
that that a n (x) = x for every x, each equivalence class has n elements. 

A further reduction can be made by considering the complement. The comple- 
ment of an n-bit binary string x = x\X2 • • • x n is x = x~\x~2 • • • x n . Note that two 
vertices x and y are adjacent iff x and y are adjacent. By considering these two 
operations, the vertices of M n is partitioned into \M n \/2n classes, each of them has 
2n vertices (Fig. [2]). Here and hereafter, we denote the number of vertices of a 



2 



p(OOOll) 




p(OOlOl) 



oo_io 
ion 



il lOOOl 



001 Oi 
1101 



p(00101) group 



Figure 2: The graph R$ and its relationship to the vertices of M5. 



graph G by \G\. 

For an n-bit binary sequence x, let p{x) denote this equivalence class including 
x, i.e., p(x) = {a 1 (x) , a 1 (x) | < i < n}. Let Rn denote the graph whose vertices 
are these equivalence classes and two vertices p{x) and p(y) in R n are adjacent iff 
there is an edge between u and v in M n for some u € p{x) and v G p{y)- 

The following lemma, which was shown by Shields and Savage [5], guarantees 
that we can lift a Hamiltonian path in R n to a Hamiltonian cycle in the middle 
levels graph. 

Lemma 1. If there is a Hamiltonian path in R n starting from the vertex p(0 +1 1 ) 
and ending at the vertex p(0(01) k ), then there is a Hamiltonian cycle in M n . 



3 Decomposition based on Runs 

Since the graph i? n is still huge (i.e., I-R37I ~ 4.8 • 10 8 ), we divide R n into a number 
of smaller graphs and search them individually and possibly in parallel. 

A run of a binary string x is a consecutive appearance of l's or 0's in x. For 
example, we say that 000000 has one run and 001011 has four runs. We will divide 
R n into three parts depending on the number of runs of strings in a vertex. Notice 
that p(x) may contain strings having different runs. We pick a string with k one's 
such that it starts with and ends with 1 as a representative of p(x), and the 
number of runs of this string is referred as the number of runs of p{x). Since this 
number is always even, we introduce a new unit called "brun" which is equal to 
two runs. 

Note that, in R n: only p(0 k+1 l k ) has 1 brun and only p(0(01) k ) has k bruns. In 
a preliminary experiment, we found that a decomposition based on the following 
three intervals is plausible (see Figs. [3] and [3]). 

• Front part : 1 ~ (|_^/2j — 1) brun(s) 
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Figure 3: The graph R 7 and the decomposition of R n based on "brun" . 



• Middle part : [k/2\ ~ (k - [k/2\ + 1) bruns 

• Rear part : (k — [k/2\ + 2) ~ k bruns 

Note that, when k = 18, these three intervals are {1, 2, . . . , 8}, {9, 10} and {11, . . . , 18}. 

We will find a Hamiltonian path in each of these three graphs and then connect 
them to get a Hamiltonian path in R n . In order to apply Lemma [TJ we fix the start 
vertex of a path in the front part to p(0 k+1 l k ) and the end vertex of a path in the 
rear part to p(0(01) fc ). In addition, we should satisfy the additional requirements 
that (i) an end vertex of a path in the front part is adjacent to a start vertex of 
a path in the middle part, and (ii) an end vertex of a path in the middle part is 
adjacent to a start vertex of a path in the rear part. 

After some considerations, we pick strings hc(k,r) := o fc ~ r + 1 (01) r T fc-r as termi- 
nals of paths. Note that p(hc(k, r)) has a maximum number of neighbors in vertices 
with r— 1 bruns and with r+1 bruns, respectively. In addition, (i) hc(k, 1) = fe+1 l fc , 
(ii)hc(k,k) = 0(01) fc , and (iii) for every i, p(hc(k,i)) and p(Rev(hc(k,i + 1))) are 
adjacent in R n where Rev(x) denotes the reverse of a string x = x\X2 ■ ■ ■ x n i.e., 
Rev(x) = x n ■ ■ ■ X2X\. We also use the following fact which can easily be verified. 

Fact 2. Let {£, £ + 1, . . . , r} be a subset of {1, 2, . . . , k}. Suppose that there is a 
Hamiltonian path in an induced subgraph of R n with vertices of at least i bruns and 
at most r bruns that starts from p(x) and ends at p(y) . Then there is a Hamiltonian 
path in the same graph that starts from p{Rev{x)) and ends at p(Rev(y)). 

For a Hamiltonian path P, let Rev(P) denote a Hamiltonian path in a same 
graph whose existence is guaranteed by Fact [2) In summary, our search procedure 



(oo 


■■ OOMOOIIOO ■■ 


■ ool 


|oo» 


■ OO] [OOOI [OO ■ 


■■OO] 



Figure 4: The decomposition of R n . Each small circle represents an induced sub- 
graph by the vertices with a specified brun. 
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is the following: First find a Hamiltonian path in each of three parts of the graph 
starting from p(hc(k,i)) and ending at p(hc(k,r)) where £ and r are the left-end 
and right-end of each interval, and let denote these three paths as Pp, Pm and Pr. 
Then connect Pp, Rcv(Pm) and Pr in this order to get a Hamiltonian path in R n 
which fulfills the condition in Lemma [TJ 

4 Ordering of Vertices 

Each vertex of the graph R n can naturally be stored using n bits of memory. 
However, this can be reduced by using an efficient ordering of the vertices. Indeed, 
since the number of vertices of i? n is less than 2 32 for n < 39, we can store them 
using a 32-bit integer par item. In this section, we give an efficient algorithm 
for ordering the vertices of our reduced graphs. A bit surprisingly, plugging this 
ordering scheme into a program gives a significant improvement of a running time 
of the program that will be shown in the next section. 

4.1 View Vertices of Middle Levels as Catalan Objects 

The n-th Catalan number is the number of expressions containing n pairs of paren- 
theses which are correctly matched and is well-known to be 



Notice that the number of vertices in R n is equal to the fe-th Catalan number C(k). 
This suggests that there is a bijection between the set of vertices of R n and the set 
of correctly matched n pairs of parentheses. 

In the following, we identify a sequence of parentheses with a binary string 
under a mapping "(" -f-> "0" and ")" o "1". In addition, by a technical reason, 
we add one "0" to the top of the string. For example, we consider that "(()(()))" 
represents the string "000100111". An 2k + 1 bit binary string starting with is 
said to be correctly matched if it is corresponding to a correctly matched n pairs of 
parentheses. 

Fact 3. For every vertex p(x) in R n , there is a unique correctly matched string in 



Proof. We should only consider a string with k one's since no string with k + 1 
one's is correctly matched. 

Suppose that we represent a string by a path in the grid such that it goes upward 
when we read and downward when we read 1. For example, a path for the string 
0000111 is drawn as Fig. [5j It is clear that a string x is correctly matched iff the 
starting point of the path for x is located at the lowest level in the path and it is 
only the point on this level. 




p(x). 
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Figure 5: A path for the string "0000111". 



Recall that p{x) contains every string that obtained from x by applying the 
cycle shift an arbitrary times. Note that, for every x with k one's, a path for x ends 
at one step higher than the starting point of the path. Hence if we draw paths for 
x and a' l (x) for some i, a path for the substring that shifted backward in a' l (x) is 
drawn at one level higher than the original level (Fig. [6|). 

cycle shift 



(a)0110100 — * 

f— (b)0001101 * 

Figure 6: A path for x = 0110100 (a) and for cr 5 (x) = 0001101 (b). A dotted line 
represents a path for the substring '01101' which goes backward by the cycle shift 
in (b). 

By this observation, it is easy to see that a correctly matched string in p(x) can 
be obtained by (i) draw a path for x, and pick the rightmost point among all points 
on the lowest level of the path, and (ii) shift x so that this point becomes the top 
of the resulting string. 

It is also easy to see that every other string in p(x) is not correctly matched. 
This guarantees the uniqueness and hence completes the proof. □ 

By this fact, there is a bijection from the set of vertices in R n to the Catalan 
objects, i.e., the vertices in R n are uniquely mapped to integers {0, 1, ... , C(k) — 1}. 

4.2 Lexicographical Ordering for Catalan Objects 

In our programs, we number vertices p(x) in R n according to the lexicographical 
ordering (starting from 0) of a correctly matched string in p(x). 

Obviously, the ordering of a string x is equal to the number of strings lexico- 
graphically smaller than x. Hence if we can count the number of strings smaller 
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Figure 7: C w (k,p) is equal to the number of left-right paths in the grid. 



than x for a given prefix x, then the ordering of x can easily be computed. For 
example, the ordering of the string 0010101 in a set S C {0, l} 7 can be computed 
as the sum of the numbers of strings in S starting from 000, 00100 and 0010100. 

Let Pi C {0, 1} 2 ^ +1 be the set of correctly matched strings of length U + 1. For 
a prefix x € {0, 1}' with t < 21 + 1, the number of strings in Pi starting with x is 
shown to be 

«■<*•'>- &T ( ?-7 ) • « 

where p = tto(^) — — 1 and k = £ — (ti(x). Here we denote the number of 0's 
and l's in x by (Jo(^) and (ti(£), respectively. Intuitively, p denotes the height of the 
end point of a path for x and k denotes the number of "remaining" one's in a string 
(see Fig. [7|). Note that these numbers are known as the Catalan Triangle (see e.g., 
the sequence A009766 of [7]). Using Eq. ([T]), we can calculate the lexicographical 
ordering of a vertex p(x) efficiently. For example, the ordering of /^(OOIOIOI) is 
given by C w (3, 2) + C w (2, 2) + C w (l, 2) = 3 + 1 + = 4. 

4.3 Runs and Narayana Numbers 

Since we decompose the graph Rn into smaller parts, it is desirable to give an 
efficient ordering algorithm for the set of vertices of these decomposed graphs. By 
a similar argument to that in Section 4.1, the number of vertices of R n with r bruns 
is shown to be 

^>-KJ)(-i)- 

which is known as the Narayana numbers. N(k, r) is the number of correctly 
matched k pairs of parentheses that contains the subsequence "()" exactly r times. 
Note that the Catalan numbers are represented by the sum of the Narayana num- 
bers, i.e., 

k 

C{k) = Y J N{k,i). 
i=l 

It is also shown that the lexicographical ordering of a string x in the set of cor- 
rectly matched strings with r bruns can be efficiently computed using the following 
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formula: 



k(k — p — r + 1) \ r J \ r — I 

that represents the number of correctly matched strings of which the meanings of 
p and k are the same as in Eq. (pQ) and r denotes the 'remaining' number of the 
subsequence "()". A detailed discussion on how to compute the ordering for such 
Narayana objects will be appeared in the full version of this note. 



5 Computational Results 

We develop a program for finding a Hamiltonian path for decomposed graphs based 
on the algorithm proposed by Shields et al.[5] in which we represent the vertices 
of graphs by the ordering described in Section UJ Using this program, we have 
succeeded to find a desired Hamiltonian path for every three parts, i.e., the front, 
middle, rear parts of R n for every 8 < k < 18, which shows the Hamiltonicity 
of the middle levels graphs for k < 18. Note that, for smaller values of k, our 
decomposition schema would not work. 

The computational results are summarized in Table [TJ Our program is executed 
on a PC with an Intel Xeon processor of 2.26 GHz and 24 GB of memory available. 
Note that the maximum memory used in our experiments was about 9 GB. We 
show the elapsed time in seconds, and the case that takes less than 1 second is 
shown as 0. 

The second column shows the elapsed time of a base program to find a path 
in the entire graph R n . In a base program, we don't use our ordering scheme and 
vertices are stored as n-bit strings. The third column shows the longest elapsed 
time of a base program for finding a path in each of three decomposed graphs. The 
fourth column shows the elapsed time of a program with the ordering technique 
for the entire graph R n . The later columns show the elapsed time of a program 
in which both techniques, i.e., the decomposition described in Section 3 and the 
ordering described in Section 4.3 are included. 

A bit surprisingly, introducing the ordering into a search program gives a sig- 
nificant improvement of the running time. The combination of our two techniques 
reduces the running time by a factor of about 30 when k = 16. For k = 18, the 
number of vertices of the front, middle and rear parts of the graph is 120, 624, 130, 
236, 390, 440 and 120, 624, 130, respectively. Notice that the running time is the 
longest for the front part of the graph. This suggests that finding a Hamiltonian 
path is harder for a graph consisting of vertices with smaller number of runs than 
that with larger number of runs. 

The source codes of the programs we used as well as some additional data are 
available on the web page [6j. Note that our program can handle up to k = 19. At 
the time of writing this note, the search for the front and rear parts of the graph 
for k = 19 has been finished successfully, and that for the middle part, which has 
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Table 1: 
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(6.3 days) 
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627,204 


511,342 


359,015 


627,204 
(7.3 days) 



about 1.18 billion nodes, is in progress. 
Note Added 

Several months after writing the above, our program for finding a Hamiltonian 
path in the middle part of the graph for k = 19 has successfully terminated. This 
confirms that the middle two levels of 39-dimensional hypercube is also Hamilto- 
nian. The number of vertices of the front, middle and rear parts of the graph is 
about 2.92 x 10 s , 1.18 x 10 9 and 2.92 x 10 s , respectively. The running time of the 
program (executed on the same machine as above) is about 56 days, 81 days and 
27 days, respectively. 
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